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We use the abstract search algorithm and its extension due to Tulsi to analyze a spatial quantum 
search algorithm that finds a marked vertex in Hanoi networks of degree 4 faster than classical 
algorithms. We also analyze the effect of using non-Groverian coins that take advantage of the small 
world structure of the Hanoi networks. We obtain the scaling of the total cost of the algorithm as a 
function of the number of vertices. We show that Tulsi's technique plays an important role to speed 
up the searching algorithm. We can improve the algorithm's efficiency by choosing a non-Groverian 
coin if we do not implement Tulsi's method. Our conclusion is based on numerical implementations. 



I. INTRODUCTION 

Grover's algorithm [l| allows one to find a marked item 
in an unsortcd database quadratically faster compared 
with the best classical algorithm. It is the paradigm for 
many quantum algorithms that use exhaustive search. 
The main technique used in Grover's algorithm, called 
amplitude amplification, can be applied in many compu- 
tational problems providing gain in time complexity. 

A related problem is to find a marked location in a spa- 
tial, physical region. Benioff [|[ asked how many steps 
are necessary for a quantum robot to find a marked ver- 
tex in a two-dimensional grid with N vertices. In his 
model, the robot can move from one vertex to an adja- 
cent one spending one time unit. Benioff showed that a 
direct application of Grover's algorithm does not provide 
improvements in the time complexity compared to a clas- 
sical robot, which is O(N). Using a different technique, 
called abstract search algorithms, Ambainis et. al. Q 
showed that it is possible to find the marked vertex with 

O (^/NlogN^ steps. Tulsi Q was able to improve this 

algorithm obtaining the time complexity O U/N log W) . 

The time needed to find a marked vertex depends on 
the spatial layout. The abstract search algorithm is a 
technique that can be applied to any regular graph. It 
is based on a modification of the standard discrete quan- 
tum walk. The coin is the Grover operator for all vertices 
except for the marked one which is —I. The choice of the 
initial condition is also essential. It must be the uniform 
superposition of all states of the computational basis of 
the coin-position space. This technique was applied with 
success on higher dimensional grids Q, honeycomb net- 
works [U , regular lattices @ and triangular networks Q ■ 
Spatial search in Hanoi network of degree 3 (HN3) was 



analyzed in Ref. [8(. Recently, the abstract search algo- 
rithm was applied for spatial search on Sierpinski gas- 
ket @. 

In this work, we analyze spatial search algorithms on 
the Hanoi network of degree 4 (HN4) extending the anal- 
ysis performed for HN3 [8| . HN4 is a special case of small 
world networks, which are being used in many contexts 
including quantum computing |10l . ITTj . We also analyze 
the use of a modified coin operator instead of the Grover 
coin, which is used in the standard form of the abstract 
search algorithm, to take advantage of the small world 
structure. Our results are based on numerical simula- 
tions, but the hierarchical structure of HN4 indicates that 
analytical results can also be obtained. Hanoi networks 
have a hierarchical structure of fractal type that helps to 
gain insights of spatial search algorithms in graphs that 
are not translational invariant. Recently, there has some 
effort in this direction [t| QjJ 

The structure of the paper is as follows. Sec. |TT] intro- 
duces the degree-4 Hanoi network. Sec. IIIII describes the 
standard coined discrete quantum walk on HN4. Sec. IIVI 
reviews the basics of the abstract search algorithm and 
Tulsi's method. Sec. |V] describes the modification on the 
coin operator we propose to enhance the time complexity 
of quantum search algorithms on HN4. Sec. I VII describes 
the main results based on numerical simulations. Finally, 
we present our final remarks in Sec. IVIII 



II. HIERARCHICAL STRUCTURES 

The Hanoi network has a cycle with N = 2™ vertices as 
a backbone structure, that is, each vertex is adjacent to 
2 neighboring vertices in this structure and extra, long- 
range edges are introduced with the goal of obtaining a 
small- world hierarchy. The labels of the vertices < k < 
2" — 1 can be factorized as 
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fc = 2 fcl (2 k 2 + l), 



(1) 
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where k\ denotes the level in the hierarchy and ki labels 
consecutive vertices within each hierarchy. In any level, 
one links the vertices with consecutive values of ki keep- 
ing the degree constraint. When fei = 0, the values of 
k are the odd integers. For HN4, we link 1 to 3, 3 to 
5, 5 to 7 and so on. The vertex with label 0, not being 
covered by Eq. (TTJ , has a loop as well as vertex of label 
2"" 1 . Fig. Q] shows all edges for HN4 when the number 
of vertices is 16. Using this figure, one can easily build 
HN4 recursively, each time doubling the number of ver- 
tices. Our analysis will be performed for a generic value 
of n to allow us to determine the computational cost as 
function of TV. 




FIG. 1: HN4 with 16 vertices. 

HN4 has a small- world structure because the diameter 
of the network only increases with ~ y/~N [IH , less fast 
than the number of vertices TV. Yet, HN4 is a regular 
graph of fixed degree d = 4 at each vertex. 



III. QUANTUM WALKS ON HIERARCHICAL 
STRUCTURES 



A coined quantum walk in HN4 with TV = 2" vertices 
has a Hilbert space H = He ® tip, where He is the 
4-dimcnsional coin subspace and Hp the TV-dimensional 
position subspace. A basis for He is the set {|a}} for 
< a < 3 and Hp is spanned by the set with, 
< k < TV — 1. We use the decomposition k — (k±, fe), 
given by Eq. (UJ), when convenient. A generic state of the 
discrete quantum walker in HN4 is 



3 N-l 



(2) 



a=0 fc=0 



The evolution operator for the standard quantum 
walk 141 is 



U = So(C®I), (3) 
where / is the identity in Hp and S is the shift operator 



defined by the following equations, 
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if k<i is even 
if &2 is odd 

if ki is even 
if ki is odd 



and 
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The arithmetical operations on the second ket is per- 
formed modulo TV. The shift operator obeys S 2 = I. C 
is a unitary coin operation in He- In the standard walk, 
C is the Grover coin, denoted by G, 



G 




(4) 



which is the most diffusive coin |15| . 

The dynamics of the standard quantum walk is given 

by 



|tf (t)> = U* |* c 



(5) 



where |^o) is th- e initial condition. After t steps of uni- 
tary evolution, we perform a position measurement which 
yields a probability distribution given by 



p k = Y^\{a,k\U t 



(0) 



o=0 



IV. ABSTRACT SEARCH ALGORITHMS 

The abstract search algorithm 0] is based on a mod- 
ified evolution operator U' = S ■ C' , obtained from the 
standard quantum walk operator U by replacing the coin 
operation C with a new unitary operation C' which is 
not restricted to He and acts differently on the searched 
vertex. The modified coin operator is 



a 



-I® |fc )(fco|+C® (I-\k )(ko\), (7) 



where kg is the marked vertex in a regular graph and C 
is the Grover coin G, the dimension of which depends on 
the degree of the graph. Ambainis et. al. 0] have shown 
that the time complexity of the spatial search algorithm 
can be obtained from the spectral decomposition of the 
evolution operator U of the unmodified quantum walk, 
which is usually simpler than that of V . 
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The initial condition \tpo) is the uniform superposi- 
tion of all states of the computational basis of the whole 
Hilbcrt space. This can be written as the tensor prod- 
uct of the uniform superposition of the computational 
basis of the coin space with the uniform superposition 
of the position space. Usually, this initial condition can 
be obtained in time 0(\/~N), where N is the number of 
vertices. 

The evolution operator is applied recursively starting 
with the initial condition \ipo). If tf is the running time 
of the algorithm, the state of the system just before mea- 
surement is U' tf \ipo)- If one analyzes the probability of 
obtaining the marked vertex k as function of time since 
the beginning of the algorithm, one gets an oscillatory 
function with the first maximum close to t / . 

After a little algebra, the evolution operator V = S-C 
can be converted into the form V = U- Rk , where Rk = 
1-2 \uc, &o) (u c , ko\, U is given by Eq. (J3J), and \u c ) is 
the uniform superposition of the computational basis of 
the coin space. Using the expression U' = U ■ Rk as a 
starting point, Tulsi proposed a new version of the search 
algorithm, which requires an extra register (an ancilla 
qubit) used as a control for the operators Rk and U. The 
operators acting on the ancilla register are described in 
Figure [2 where — Z is the negative of Pauli's Z operator 
and 



X s = 



cos 6 sin S 
— sin 5 cos 5 



(8) 



The value of S is the one that optimize the cost of the 
algorithm. For two-dimensional lattices, Tulsi Q showed 
that cosS oc 1/Vlog N. 
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FIG. 2: Tulsi's circuit diagram for the one-step evolution 
operator of the quantum walk search algorithm. \up) is the 
uniform superposition of the computational basis of the 
position space. 

The Tulsi's evolution operator is 



U" 



{-Z®I)- C{U) ■ {x\ ® I) ■ C(R kQ ) ■ (X s ® I), (9) 



where C(U) and C(Rk ) are the controlled operations 
shown in Figure[2]and / is the identity operator in T-L. We 
want to determine how many times U" must be iterated, 
taking |1) \uq) \up) as the initial condition, in order to 
maximize the overlap with the search element. 



V. MODIFIED METHOD 

The coin in a quantum walk is used to determine the di- 
rection of the movement. The Grover coin is an isotropic 
operator regarding all outgoing edges from a vertex. It is 
useful in networks that have no special directions, such 
as two-dimensional grids and hypercubes. The Hanoi 
network, on the other hand, has a special direction that 
creates the small world structure. Any edge that takes 
the walker outside the circular backbone provides an in- 
teresting opportunity in terms of searching. The strategy 
is to have a parameter that can control the probability 
flux among the edges, reinforcing or decreasing the flux 
outwards or inwards the circular backbone. 

Instead of using the Grover coin of the abstract search 
algorithms, we analyze the use of a modified coin given 
by 

C=§(|0) (0| + |1)(1| + |0)<1| + |1) (0|) + 



I e{d - 2e) 
<P{d~2) 



d-l 



Z>'> (o| + |o) 01) 



J=2 
1 



(10) 



J=2 



d 2 (d-2) 

2(d-2e) 
d(d-2) 



d-i 



E (li>0"l)-J. 



where d = 4 is the degree at each vertex. When e = 1, 
the Grover coin is recovered. When < e < 1, the 
probability flux along small- world edges (labels and 1) 
which escapes from the circular backbone is weakened. 
When 1 < e < 4, the probability flux off the backbone 
is reinforced, allowing the walker to use the small world 
structure with higher efficiency. Hence, this new coin 
controls the bias to escape off the circular backbone of 
HN4 through the parameters e. 

The abstract search algorithms use a uniform distri- 
bution as initial condition. We change this recipe. The 
initial condition is 




l#>)> = \/ 5 (|0>l*> + |l>l*>)- 



3) \s) , 
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d(d-2)£ 2 



where \s) is the uniform superposition on the position 
space. When e = 1 the initial condition is the uniform 
superposition of coin-position space. 

We want to check whether it is possible to improve the 
abstract search algorithm by modifying the coin operator 
for HN4 in such way we can tune parameter e for obtain- 
ing the best rate of probability flux between the circle 
backbone and small- world edges. In a previous paper [8|, 
we have concluded that, for HN3 without using Tulsi's 
method, it is better to choose e = 1.75. After analyzing 
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further this issue, we have to reconsider this conclusion, 
mainly when one uses Tulsi's method, which seems to 
favor the Grover coin even in nonhomogeneous graphs. 

In this work, we consider three different methods: 1) 
the abstract search algorithm, 2) the Tulsi's method, and 
3) the modified method. The analysis of the evolution of 
the quantum search algorithm using the new coin and ini- 
tial condition is far more complex than the standard one. 
Our conclusions here are based in numerical simulations. 



VI. MAIN RESULTS 



method (lower curves) and the abstract search algorithm 
using Tulsi's method (higher curves). The curves are 
very flat around e = 1, which correspond to the Grover 
coin. This shows the the modified method does not play 
an important role for improving the efficiency of the al- 
gorithm. These curves do not provide enough clues for 
choosing e. The final answer can be achieved by ana- 
lyzing the effect of e on the total cost of the algorithm. 
We measure the cost as the number of times the evolu- 
tion operator is applied, or equivalently the number of 
oracle queries considering the repetitions necessary for 
amplitude amplification. 



Fig. [3] shows the oscillatory behavior of the probabil- 
ity of finding the walker at the marked vertex using the 
modified method with e = 0.75 (lower curve) and the 
abstract search algorithm using Tulsi's method (higher 
curve). Initially, the probability is close to zero, because 
the initial condition is a state that is close to the uni- 
form superposition of all vertices. The running time of 
the algorithm is the value of t for which the probability 
is close to its first maximum. Note that, without using 
Tulsi's method, the maximum value of the probability is 
smaller than that of the abstract search algorithm with 
Tulsi's method. In either case, the maximum value of 
the probability is not close to 1, as one would expect in 
order to have high probability to find the marked vertex. 
This means that the algorithm must be rerun many times 
to amplify the success probability. For the lower curve, 
the number of repetitions is large, in fact, it scales with 
N, which has a strong impact on the total cost of the 
algorithm. We call success probability the value of the 
probability for the first peak of Fig. [3J The marked ver- 
tex used in all simulations is k = 3, but the conclusions 
will not depend on the hierarchy of the target vertex. 




20 40 60 80 100 120 140 160 180 200 
steps 

FIG. 3: Probability of finding walker at vertex k = 3 as a 
function of time using the modified method and the abstract 
search algorithm with Tulsi's method. 

Now let us try to answer the following question about 
parameter e: What is the best value of e for the spatial 
search algorithm? Fig. @] shows the success probability as 
function of e for three values of N both for the modified 
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Tulsi method, N=2048 
Tulsi method, N=4096 
modified method, N=1024 
modified method, N=2048 
modified method, N=4096 
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FIG. 4: Success probability of the search algorithm as 
function of e for three values of N using both the modified 
coin and Tulsi's method on the top of the abstract search 
algorithm. 

Fig. [5] shows the total cost of the search algorithm as 
function of e for two values of N both for the modified 
method (higher curves) and the abstract search algorithm 
using Tulsi's method (lower curves). The curves are very 
flat around e = 1 as before, but we can conclude that the 
best values are e = 0.75 for the modified method and e = 
1 (Grover coin) for the abstract search algorithm using 
Tulsi's method. From now on, we will take these values 
of e for the rest of this paper. The fact that the algorithm 
cannot be improved by modifying the coin after Tulsi's 
method seems to show that the algorithm has achieved its 
best performance using the Grover coin. This conclusion 
for a graph that has non-homegeneous vertices and non- 
isotropic edges was not the one expected by us at the 
beginning and it is quite surprising. 

Fig. [5] shows the success probability after a single run 
of the search algorithm as a function of the network size 
N in log-scale for the modified method (lower points) 
and the abstract search algorithm using Tulsi's method 
(higher points) . From the inclination of the best fitting 
line, we conclude that the success probability of the mod- 
ified method decays approximately as 0.62/iV ' 37 . Using 
the technique of amplitude amplification [H, the algo- 
rithm must be rerun around O(iV 0185 ) times in order to 
ensure a final probability close to 1. This produces a high 
impact on the total cost of the algorithm. Recall that for 
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overhead of 1/y/P, where P = 0{l/N° 37 ) (sec Fig. [BJ. 



FIG. 5: Computational cost of the search algorithm in terms 
of e for two values of N using both the modified coin and 
Tulsi's method on the top the abstract search algorithm. 



the two-dimensional grid, the number of repetitions is 
O(yiogTV) 0. This result shows that Tulsi's method 
plays an important role in terms of computational com- 
plexity. The best value of 5 in Eq. {5} for HN4 seems 
to be cos S = 0(1/ log N). In Fig. [6j we see that success 
probability after a single run of the abstract search al- 
gorithm using Tulsi's method does not depend on N. In 
this case, we can rerun the algorithm a fixed number of 
times to obtain an overall probability very close to 1. 
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FIG. 6: Success probability of the search algorithm as a 
function of N. The fitting shows that P = 0.62/iV - 37 for 
the modified method. 

Fig. [7] shows the computational cost of the search al- 
gorithm as a function of the network size N both for the 
modified method (cross points) and the abstract search 
algorithm using Tulsi's method (x points). We have not 
used the method of amplitude amplification in this exper- 
iment. We have displayed the best fitting lines for both 
cases, which scale as O(N 65 ). For the Tulsi method, 
this is the total cost of the algorithm, because the suc- 
cess probability is high (see Fig. [5]). For the modified 
method, the total cost is O(iV - 84 ), because we must use 
the method of amplitude amplification, which puts an 
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FIG. 7: Computational cost of the search algorithm in terms 
of N in log scale. The dotted curves are fitting curves. For 
the modified method, the best fitting is 1.25iV ' 65 . For the 
Tulsi method, the best fitting is 2.62JV ' 65 . 

It is important to note that the data in Fig. [7J is not 
precise enough to detect the presence of log N terms in 
the expression of the total cost. We are using results 
with small N to draw conclusions about the asymptotic 
behavior. There are imprecisions in the simulations that 
come from the discrete nature of the spatial layout. For 
instance, the lower curve in Fig. [3] has quick oscillations 
that have some impact on the impreciseness of the total 
cost. Usually, the results using Tulsi's method are more 
stable. 

The scaling of the total cost has a slight variation when 
we change the position of the target. On the other hand, 
the prefactor changes. In terms of graph structure, HN4 
is non-homogeneous, because the vertices can be divided 
in hierarchical levels. This non-homogeneity does not 
play an important role in term of the cost of finding a 
vertex. The same kind of conclusions holds for HN3, 
analyzed in Ref. Q. 

We have redone and extended the simulations for HN3, 
using a new implementation. The scaling for the to- 
tal cost in terms of N is (1) O(N 62 ) using Tulsi's 
method, and (2) O(N - 74 ) using the modified method 
with e = 1.75 applying amplitude amplification. There 
are two important conclusions we draw from the com- 
parison between HN3 and HN4, one is about the graph 
degree and the other about the value of e. 

The degree of the graph seems to play no important 
role in terms of efficiency Similar conclusions were draw 
in Ref. 0, which compared the efficiency of quantum 
search algorithms in triangular, square, and hexagonal 
lattices that have degrees 3, 4, and 6, respectively. The 
scaling of the cost as a function of TV is the same for all 
of them. 

The optimal value of e for HN3 is larger than 1 using 
the modified method, which means that the probability 
flux toward the edges leaving the circle backbone is en- 
hanced. For HN4, the optimal value of e is 0.75, smaller 
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than 1. We cannot say that long range connections (in 
term of hierarchical level) helps in the quantum search. 
In fact, for HN4 we have to decrease the probability flux 
in the edges leaving the circle backbone. This can be in- 
terpreted, when we take into account that HN4 is really 
small-world and mean-field like in terms of the average 
distance between any two vertices, which scales logarith- 
mically with system size, whereas the average distance 
scales as y/~N for HN3. That means, in HN4 it is less sig- 
nificant to take long-range jumps, because a random mix 
is already enough to get to most other sites; whereas in 
HN3, if the walker does not take more long-range jumps 
than nearest-neighbor jumps, it is difficult to go very far. 



VII. FINAL REMARKS 

We have analyzed spatial search algorithms on degree- 
4 Hanoi networks with the goal of extending the abstract- 
search-algorithm technique for nonhomogencous graph 
structures with fractal nature. We have proposed a mod- 
ification of the abstract search algorithm by choosing a 



coin that takes advantage of the edge asymmetry of HN4. 
We have obtained a faster algorithm by tuning numeri- 
cally parameter e. The cost of this algorithm is O(N 84 ) 
in terms of number of oracle queries. The algorithm uses 
the standard method of amplitude amplification on top 
of the modified abstract search algorithm with e = 0.75. 
This value of e tells us that the probability flux is higher 
on the circle backbone of HN4 than on the edges that 
produces the small world structure. 

We have also analyzed Tulsi's method on top of the 
abstract search algorithm. In this case, the Grovcr coin 
(e = I) seems to be the best option and the modi- 
fied method does not improve the algorithm. The cost 
of the algorithm is O(N 65 ). This is above the lower 
bound, which is O(N 05 ), and above the cost of search- 
ing a marked vertex on two-dimensional lattices, which 
is 0((N logiV) - 5 ). We have used numerical methods to 
estimate the cost scale. This means that logiV factors 
may be lost, which could decrease the scale of 0.65 in the 
total cost. 

Our works on progress are now focused on obtaining 
analytical results regarding search algorithms and general 
quantum walks for the Hanoi network HN4. 
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